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1 Introduction 



The canonical height h on an abeUan variety A defined over a global field 
k is an object of fundamental importance in the study of the arithmetic of 
A. For many applications it is required to compute h{P) for a given point 
P € A{k). For instance, given generators of a subgroup of the Mordell-Weil 
group A(k) of finite index, this is necessary for most known approaches to the 
computation of generators of the Mordell-Weil group A{k). Furthermore, the 
regulator of A{k), which appears in the statement of the conjecture of Birch and 
Swinnerton-Dyer, is defined in terms of the canonical height and thus we need 
the ability to compute canonical heights in order to gather numerical evidence 
for the conjecture in the case of positive rank. 

Here we are concerned with the case where A is the Jacobian variety of a smooth 
projective curve C of genus g over k. If f/ < 3, it is known how to compute 
canonical heights using arithmetic on an explicit embedding of the Kummer 
variety K of A into P^""! (cf. [3i [HI 1371 [MI and [gSl Chapter 3]). Trying 
to imitate this approach in the higher genus case quickly causes problems, as 
the Kummer variety becomes rather complicated (see the discussion in |28l 
Chapter 4]). 

Instead we propose to use a result due to Faltings [13] and Hriljac [H] (see The- 
orem 13. 2p , expressing the canonical height in terms of arithmetic intersection 
theory. Here the non-archimedean intersection multiplicities take place on regu- 
lar models of C, whereas the archimedean intersection multiplicities are given in 
terms of Green's functions on the Riemann surface associated to C. In Section 
[5] we discuss the local theory before putting the local results together in Section 
[31 culminating in Theorem 13.21 which establishes the connection between h and 
arithmetic intersection theory. 

In Section[4]we show how the necessary arithmetic intersection multiplicities can 
be computed in practice. In the non-archimedean case we reduce the problem to 
the computation of certain Grobner bases. Then we show that the archimedean 
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intersection multiplicities can be computed using theta functions with respect to 
the complex torus C^/A associated to A. In order to make these steps practical, 
we need to be able to decompose divisors into prime divisors and to work on 
C^/A explicitly. 

We present a practical algorithm, implemented in the computer algebra system 
Magma [25| , for the computation of h for hyperelliptic curves in Section [5] by 
explaining how these two points can be resolved in that case. Several examples 
are given in Section |51 where the performance of the algorithm is investigated 
as well. Finally, we elaborate on what is needed to extend the algorithm to non- 
hyperelliptic curves. A different, but similar, algorithm for the computation 
of h using arithmetic intersection theory has been developed independently by 
Holmes [20 . 

Acknowledgments: The research presented here is also described in Chapters 5 
and 6 of my PhD dissertation [28]. I would like to thank my advisor Michael 
StoU for suggesting this topic and for his constant help and encouragement. 

Some of the research described in this work was conducted while I was visiting 
the University of Warwick and the University of Sydney. It is my pleasure to 
thank both institutions for their hospitality, as well as David Holmes, Samir 
Siksek and Steve Donnelly for their invitations. 

2 Local Neron symbols 



In this section we discuss the theory of local Neron symbols whose existence 
was first proved by Neron in [30) . We shall present an interpretation that is 
suitable for explicit computations, following essentially Gross [T7] and Hriljac 
[21] . The content of the latter work is also discussed by Lang in |23l . In order to 
present these results, we need to briefly recall some basic notions of intersection 
theory on arithmetic surfaces. 

In the following 3 sections C denotes a smooth projective geometrically con- 
nected curve of positive genus g, defined over a field k which will be specified 
as we go along. Let Div(C) denote the group of divisors on C Xj, k^°P, where 
k^°P is a separable closure of k. For an extension k' of k contained in k^°P we 
denote the subgroup of fc'-rational divisors by Div(C)(fc'). For each n G Z the 
set Div"(C) is defined to be the set of divisors of degree equal to n and we set 

Div"(C)(fc) Div"(C) n Div(C)(fc). 

If / e fc(C)* and D ^ Ejm-jiQj) £ Div"(C)(fc) is relatively prime to div(/), 
then we set f{D) := JJ. f{QjT'. 

Let /c be a non-archimedean local field valuation v with discrete valuation ring 
R, uniformizing element tt, residue field £ and spectrum S — Spec(i?). 

Definition 2.1. A model ip : C ^ S oi C over S is an integral, flat, Noetherian 
S'-scheme of dimension 2 whose generic fiber is isomorphic to C. 

Let ijj : C ^ S denote a model of C over S. By abuse of notation, we often omit 
ip and simply call C a model of C over S. We denote the special fiber of C by 
Cv Then is connected by [31| Corollary 8.3.6]. 
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Let Div(C) denote the group of Weil divisors on C. If D £ Div(C)(fc) is prime, 
then we write Dq for the Zariski closure of -D on C. This is a prime divisor on 
C and we extend the operation D i— > Dc to all of Div(C)(A:) by linearity. 

We want to use intersection theory on models of C over S. Although this can 
be defined more generally, it is convenient to restrict to proper regular models. 
For a proof that such a model always exists, see [2H §8.3.4]. In our algorithm 
we shall need a proper regular model of a specific kind; this will be discussed in 
Section 1431 

So suppose that ip : C — > 5 is a proper regular model of C over S. 

Definition 2.2. [231 §111.2] Let D, E be two effective divisors on C without 
common component and let P G C^, be a closed point. Let Id,p and Ie,p be 
defining ideals of D and E, respectively, in the local ring Oq.p- Then the integer 

ip{D,E) lengtho^^ {Oc,p / {Id^p + Ie,p)) 

is called the intersection multiplicity of D and E at P. The total intersection 
multiplicity of D and E is 

i,{D,E) ■.= Y,zp{D,EmP):l], 

p 

where the sum is over all closed points P d Cy. Finally, we extend i p and iy by 
linearity to divisors D,E d Div(C) without common component. 

A fibral Q-divisor is an element of the Q- vector space Divt,(C) generated by the 
irreducible components of C^,. If D G Div°(C)(fc), then we denote by ^v.c{D) 
a fibral Q-divisor on C such that Dc + ^v^c{D) has trivial intersection with all 
elements of Div„(C). That such a fibral Q-divisor always exists was first proved 
by Hriljac, cf. ^ Theorem III.3.6]. 

Now we have assembled all ingredients necessary to define the central objects 
of this section in the non-archimedean case. 

Definition 2.3. The local Neron symbol on C over k is defined on divisors 
D,E £ Div"(C)(fc) with disjoint support by 

{D,E), ■.^i,iDc + <i>y^c{D),Ec)\og#t 

Remark 2.4. [221 Theorem III. 5. 2]. The local Neron symbol depends neither on 
the choice of the regular model C nor on the choice of ^^^ciD). 

Next we let k denote an archimedean local field. We can assume k = C (see 
Proposition I2.8f d]) below), so that C{k) is actually a compact Riemann sur- 
face. In arithmetic intersection theory one uses Green's functions to define 
archimedean intersection multiplicities, but for us a somewhat weaker notion 
suffices. The next result follows from [251 Theorem 13.5.2] combined with [231 
Proposition II. 1.3], see also [17, §3]. 

Proposition 2.5. Let X be a compact Riemann surface and let dji be a positive 
volume form on X such that J^d^j. = 1. For each E G Div(X) there exists a 
function 

gs ■■ X\sviY>T>{E) -> M, 

called an almost-Green's function with respect to E and d^, such that the fol- 
lowing properties are satisfied: 
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(i) The function qe is C°° outside o/supp(£') and has a logarithmic singu- 
larity along E. 

1 

deg{E)d^jL = —ddgE- 

TT 

Let V be the absolute value on k and fix a volume form rf/i on C(fc) such that 
Ix = 1- 

Definition 2.6. The pairing {■,-)v that associates to all D,E E Div''(C)(fc) 
with disjoint support the intersection multiplicity 

iy{D,E) -.^gEiD) 

is called the local Neron symbol on C over k. 

Remark 2.7. It follows from [221 Proposition II. 1.3] that the local Neron symbol 
does not depend on the choice of gE or dfi. See also [231 Theorem III. 5. 3]. 

We list the most important properties of the local Neron symbol, both non- 
archimedean and archimedean, in the following proposition. 

Proposition 2.8. (Neron, Gross, Hriljac) Let k be a local field with valua- 
tion V. The local Neron symbol satisfies the following properties, where D,E E 
Div''(C)(fc) have disjoint support. 

(a) The symbol is bilinear. 

(b) The symbol is symmetric. 

(c) If f e k{C)*, then we have {D,div{f))^ ^ -log \f{D)\y. 

(d) If k' is a finite extension of k with valuation v' extending v, then we have 
{D,E)„, = [k' : k]{D,E),. 

Proof See [H §111.5] and [21 Theorems 11.3.6, 11.3.7]. 

□ 

3 Neron symbols and canonical heights 



In this section we let k denote a global field with ring of integers Ok . We assume 
that C is given by C^-integral equations. Let Mk = M^UM^ denote the set of 
places of fc, with absolute values | ■ |i, normalized to satisfy the product formula. 
Here (respectively M^) denotes the set of non- archimedean (respectively 
archimedean) places. For each place v £ Mk we let k^ denote the completion of 
k a,t v. If w G M°, we let Oy be the ring of integers at v. 

Let A denote the Jacobian variety of C and let K denote its Kummer variety 
A/{±1}. Let K ^ P^""! be an embedding of K and let 

K : A — » K ^ p2'-i. 
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Definition 3.1. [SO] , [T7] , [TS] The canonical height (or Ncron-Tate height) on 
A is the function defined by 

h{P) lim 4-"/i(k(2"P)), 

n— J-oo 

where h is the usual absolute height on P^'^^. The canonical height pairing (or 
Neron-Tate height pairing) on A is defined by 

{P, Q)nt \{hiP + 0) - h{P) ~ h{Q j). 

Note that taking the absolute height on P-^"^^ means that h does not depend 
on k. 

Next we shall relate the canonical height to Neron symbols. If £) G Div(C)(fc) 
and V e Mk, then we define Dy := D (»k K- For D,E & Div°(C)(fc) with 
disjoint support and v G Mk we define {D,E)y :— {Dy,Ey)y and the global 
Neron symbol by 

(D,E) (^'^)- 
veMk 

This is a finite sum, since over all places of good reduction the Zariski closure 'C" 
of the given equations of C fc„ over Spec(Ot,) is a proper regular model over 
Spec(C'^). Hence we have $ J, (D„) = for all such w and ^^(Z)^ , P") ^0 
for only finitely many such v. 

By Proposition I2.8f c) and the product formula, {D, E) only depends on the 
classes [I?], [E] G Pic°(C) and hence we can drop the assumption that D and E 
have disjoint support. 

Tiieorem 3.2. (Faltings Hriljac \21f ) Suppose C is a smooth projective 
geometrically connected curve of positive genus g defined over a global field k. 
IfD,Ee Div°(C)(A:), then we have 

(D,E) = -{[D],[E])^T. 

The practical importance of this result lies in the fact that we can, at least in 
principle, compute the canonical height on the Jacobian using data associated 
to the curve. We do not impose any further conditions on C (yet). Suppose 
that we are given a point P G A{k) and we want to compute its canonical height 
h{P). In order to use Theorem 13.21 for this purpose, we proceed as follows: 

(1) Find divisors D,E e Div°(C)(fc) such that [D] ^ [E] ^ P and supp(D) n 
supp(£') = 0. 

(2) Determine a finite set U of places v G such that {v G M^^ : {D,E)y ^ 
0} C U. 

(3) Find a proper regular model C oi C (E)k ky over Spec(O^) for all v E U. 

(4) Compute iv{Dv,CiEv.c) for all v eU. 

(5) Compute iv{^v.c{Dvfi)i Ev,c) for all i> G J7 of bad reduction. We call this 
the correction term. 
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(6) Find an almost-Green's function g^^, and compute gE^{Dy) for all v G 

(7) Sum up all local Neron symbols. 

We deal with these steps in the following sections. 

Remark 3.3. We shall tacitly assume from now on that step (1) is always possible 
in principle, that is every P we encounter can be represented using a /c-rational 
divisor. According to [331 Proposition 3.3] this is guaranteed whenever the curve 
has a fc„ -rational divisor of degree 1 for all v G Mk- If we have P G A{k) which 
cannot be represented using a fc-rational divisor, then we have two options: 

1) Work over a field extension k' of k such that there exists some D £ 
Div"(C)(fc') satisfying [D] = P. 

2) Compute a multiple nP such that there exists D G Div*'(C)(fc) satisfying 
[D] = nP and use the quadraticity of the canonical height. 

The existence of n as in 2) follows from [33l Proposition 3.2]; we can take for n 
the period of C over k. 

4 Computing Neron symbols 



In this section we shall address the steps needed for the computation of global 
Neron symbols introduced in the previous section. The first two steps are global 
in nature and can be viewed as preparatory steps for the remaining four sections 
which are local. 

4.1 Finding suitable divisors of degree zero 

The basic reference for large parts of the remainder of this section is 19 . If an 
ideal / is generated by elements 6i, . . . , 6„, then we write / = . . . , 6„). Let 
k be an arbitrary field. There are essentially two ways to represent a divisor 
D e Div(C)(fc). 

(a) As a sum 

i 

where Di G Div(C)(fc) is irreducible over k and to,; G Z for all i. We call 
this the free representation of D. 

(b) Assuming D is effective, using a defining ideal 

lDCk[C% 

where C° is an affine chart of C containing D. We call this an ideal 
representation of D. 
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Since in our intended applications we are allowed (and often even required) to 
vary divisors in their linear equivalence classes, it is a natural question to ask 
whether it is possible to find divisors linearly equivalent to a given divisor in a 
way that facilitates explicit computations. 

Lemma 4.1. (Hess) For all D E Div(C)(fc) and effective A e Div(C)(fc) there 
exists an effectively computable triple {D, r, a), where D G Div(C)(/c) is effective, 
r G Z and a G k{C) such that deg{D) < g + deg{A) and we have 

D = D + rA + div(a). 

We call D a reduction of D along A. // deg(A) = 1, then D is the unique 
effective divisor such that dini(£(_D — r' A)) = for all r' > 1. In this case we 
have D ^ D + rA, where r E 1 is the maximal integer such that dini(£(_D — 
rA)) = 1. 

Proof See [HI §8]. □ 

Now assume that fc is a global field, that we are given some divisor D G 
Div"(C)(fc) and we want to find E ^ D such that E and D have disjoint 
support. In other words, we are looking for an effective version of the moving 
lemma. However, we would like to keep the computations as simple as possible 
and this means that we would like to work with divisors that are reduced along 
some effective divisor of small degree whenever possible. 

This leads to the following method: 

1. Pick two effective divisors A, A' G Div(C)(fc) with disjoint support. 

2. Compute multiples nZ?, where n = 1, —I, 2, —2, . . . and reduce them along 
A and A' until we find some n and n' such that the reduction Z)„ of nD 
along A and the reduction D„i of n'D along A' have disjoint support. 

3. Let rn,rn' G Z such that nD ^ Dn + rnA and n'D ^ D'n' + rn'A' . 
Compute 

{D,D) = -^{Drr + r^A,D\,, +r.n.A') 
nn' 

^ -1^{D^,D'„,) + ^{A,D\,) + ^{D^,A') + ^-^{A,A'). 
nn' nn' nn' nn 

In practice integers n, n' of fairly small absolute value usually suffice. 

Remark 4.2. In the method above, it is not obvious how to pick A and A' in a 
way that facilitates explicit computations. If we have fc-rational divisors A, A' 
of degree 1 on C then they can be used. If C is a plane curve, then we can 
use the zero or pole divisors of functions of the form x — C, where C G fc. See 
Section [Ql for the case of hyperelliptic curves. In general the choice of A and 
of A' depends on the specific situation. 



7 



4.2 Determining relevant non-archimedean places 



We continue to let k denote a global field. Given two divisors D and E with 
disjoint support, we have to find a finite set of places v G such that any 
non-archimedean v satisfying {D,E)y ^ must lie in U. 

We can assume that D and E are effective and use their respective ideal rep- 
resentations. The idea is to cover our curve by affine patches , . . . , and 
determine the relevant places for each patch using Grobner bases. We refer to 
[TJ Chapter 4] for an introduction to the theory and applications of Grobner 
bases for polynomial rings over Euclidean rings. 

So let 

C* = Spec k[xi, . . . , Xn]/{Gi^i{xi, . . . , Xn), ■ ■ ■ , Gi.rrii {xi, ■ ■ ■ , Xn)) 

be such an affine patch, where Gij(xi, . . . , x„) S Ok[xi, . . . , Xn] for all j. From 
now on we will assume that for each v € M° there is some i,j such that Gij 
is a w-adic unit. Note that this implies that the Zariski closure Cj" of C Xj, fc„ 
over Spec(C'„) is a model for C over Spec(C'u). 

Suppose for now that the ring of integers Ok is Euclidean and that D and E 
are represented by ideals lD,i and lE,i, respectively, on C* for each i. In fact 
we can assume that lD,i and lE,i are given by bases whose elements are in 
Ok[xi, . . . , Xn]- If we compute a Grobner basis Bi of 

lD,E,i ■= (Gi,i(a;i, . . . , Xn), • ■ • , Gi^„i-{xi, . . . , X„)) -t- lD,i + lE,i 

over Ok, then Bi contains a unique element qd.ea G Ok- 
We define the set U by 

U ■-— {v e ■- OTdv{qD.E.i) > for some i}- 

For the proof of the following Lemma, we need the notion and existence of a 
desingularization in the strong sense of C for each v- The necessary details are 
presented in Section H31 below. 

Lemma 4.3. Any non-archimedean place v such that {D,E)^ > is contained 
in U . 

Proof- Suppose that {D, E)y > and let ^ : C — > C" denote a desingulariza- 
tion of C in the strong sense. The existence of ^ is asserted by 14.51 since by 
assumption C is a model of C Xk ky over Spec(C'i,). Then we must have 

iv{^v,c{D),E,^c) > and iv{Dyfi,^v,c{E)) > (1) 

or 

iviD,^c,E,^c)>0- (2) 

If (HI) holds, then v must be a place of bad reduction such that both -^j^ 
and E^yj^ intersect the singular locus of , since otherwise either <^„fi(D) or 
'^vfiiE) vanish. 
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If dD) holds, the fact that ^ is an isomorphism outside the singular locus of C 
implies that the closures -^^ and do not have disjoint supports. But 

this means that there is a point in the support of D and a point in the support of 
E having the same reduction modulo v. The claim follows easily from this. □ 

Hence the problem of determining U comes down to a combination of computing 
Grobner bases and factoring. 

If Ok is not a Euclidean ring, then we can still use this Grobner basis approach 
by writing k as fc'(a) for a primitive element a of fc over k', where fc' = Q if 
fc is a number field and k' = Fp(T) if char(/c) = p ^ 0. This trick appears in 
[TJ Exercise 4.3.1]. We add a new variable t to Ok'[xi, . . . ,a;„], satisfying the 
relation 

where (pa is the minimal polynomial of a over k' , and replace any occurrence 
of a in lD,E,i by t. Now we get at most one qD,E,i{t) G Ok'[t] \ Ok' in the 
Grobner basis of lD,E,i, but we might also have some q'u e i ^ ^k' ■ We factor 
the principal ideal {qD,E,i{ct)) in Ok and, if necessary, the principal ideal (g^ ^ ^) 
in Ok to find the relevant v G M". 

4.3 Computing regular models 



In the following three sections we let k denote a non-archimedean local field with 
valuation v. Let R be its discrete valuation ring with spectrum S = Spec(i?), 
uniformizing element tt and residue field t. Suppose that C is given by i?-integral 
equations and that the Zariski closure C of C over S' is a model of C over S. In 
practice, it is easy to find such equations for C (for instance by requiring that 
there is some unit amongst their coefficients). 

Definition 4.4. Let C be a model of C over S. A desingularization ofC in the 
strong sense is a proper birational morphism ^ : C — C such that C is regular 
and ^ is an isomorphism above every regular point of C. 

The proof of the following result can be found in Corollary 8.3.51]. It 
extends a proof due to Lipman, recalled in [2] , of the same result in the special 
case that the given model of C over S is excellent. 

Lemma 4.5. LetC he a model of C over S. Then there exists a desingulariza- 
tion of C in the strong sense. 

Although the theory works for any proper regular model of C over S, for our 
algorithm we need a desingularization C ^ C in the strong sense, as in the 
previous section. 

Lipman's proof is effective: The idea is to normalize C and then blow up the 
resulting model along its (necessarily isolated) irregular points. Repeating this 
process yields a desingularization of C in the strong sense after finitely many 
steps. The main problem is that normalizations are more difficult from a com- 
putational point of view than blow-ups. 
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A different effective proof of the existence of a desingularization of C in the 
strong sense is given by Cossart, Jannsen and Saito in [7], although they assume 
that C is excellent. The advantage of their method is that it only involves blow- 
ups, either along isolated irregular points or along irreducible components if the 
singular locus of the respective model has positive dimension. 

This method has been implemented in Magma by Donnelly. The data that can be 
accessed once C has been constructed using Magma includes the blow-up maps on 
enough afhnc patches to cover all intermediate models, the intersection matrix 
of Cy and the multiplicities of the irreducible components. 

A subtle point is that in the proof in [7 the blow-ups have to be performed in a 
specific order and this has not yet been included in the implementation, but it 
should be possible without too much difSculty. In any case, the implementation 
works well in practice. The only essential restrictions at the moment are that the 
curve has to be planar and that blow-ups along components are not implemented 
unless C is defined over Q. 

4.4 Computing non-archimedean intersection multiplici- 
ties 



We keep the notation from the previous section and fix a desingularization 
^ : C C in the strong sense, covered by affine patches 

C = Speci?[a;i, . . . ,Xs^]/ {Hi^i{xi, . . . . . . , iJ^,*, (xi, . . . ,a;sj). 

For computational purposes we shall assume for the moment that we have two 
effective divisors D and E with disjoint support whose closures Dc and Ec lie 
entirely in an affine piece C*. 

The following lemma is a well-known result from commutative algebra saying 
that quotients and localizations commute. 

Lemma 4.6. Let A be a commutative ring with unity and let T G A be a 
multiplicative subset. Let I G A be an ideal and let T denote the image of T in 
A/L. Then we have 

At I L At = {AlLA)f, 
where the subscripts denote localizations. 

Proof. See Theorem 4.2]. □ 

Let lD,i and lE,i denote defining ideals of Dq and Ec in the ring Ogi, respec- 
tively. For the computation of the intersection multiplicity we use the following 
version of the Chinese remainder theorem for modules. 

Proposition 4.7. Let A be a commutative ring and let M be an Artinian and 
Noetherian A-module. Then there is an isomorphism of A-modules 

p 

where the sum is over all maximal ideals P of A and Mp denotes the localization 
of M at P. 
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Proof. See [H Theorem 2.13]. □ 

Proposition 4.8. Suppose that Dq H Ec only intersects a single component T 
of CI- Let Or denote the ring of regular functions on T. Then we have 

iv{Dc,Ec) = lengthen {0^/{ID^^ + IE,^)0^) 

Proof. From Proposition 14. 71 we get an isomorpliism of Or-niodules 

0^/{ID,^ + lEa) = Oc%p/(^D,. + Isa), (3) 
P 

wliere tlie sum is over all maximal ideals of Or, that is, over all closed points 
P G r. By our assumptions we have 

i^iDcEc) -^^p(I?c,^^c)[^(^) :«] 
p 

= 5]lengtho^^^ {Oc.,p/{Id.+Ie.)) [KP) : «] 
p 

= lengthen (Oc.,p/(/d,. + IE,^)) 
p 

= lengtho, (Ocvp/(/d,. + lE,i)) 
p 

= lengthen (Or/(/i5,, 

using (jSj, additivity of the length and the fact that if Ad is an Op-module that 
is also an O^i p-module for some closed point P £ F, then we have 

lengtho,(M) = lengtho^^^ (M)[«(P) : I]. 

□ 

Instead of computing length^,^ {Or/{lD,i + Iea)Oy-) for each component F of 
C^, we can proceed more directly. Let 

AD,E,i,v ■= iR[xi, . . ■ ,Xs,]/lD.E.i,v)(^^) (4) 

where 

Id,E,Lv = {Hi{xi, . . .,Xsi) . ■ .,Hti{xi,. . .,Xs,)) + lD,i + lE,i- (5) 

Corollary 4.9. We have 

iy{Dc,Ec) = lengthg^^ Ad.e.z.v 

Proof. Use Proposition 14.81 and additivity of the length. □ 



Computing length^^. Ap)^E.i,v is rather easy and can be done, for instance, in 
Magma. The crucial step is the computation of a Grobner basis B of Id, E,i,v over 
the Euclidean ring R. Here the property of B that we need is that for every 
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h G R[xi, . . . ,Xsi] multivariate division of /i by -B yields a unique remainder 
that we call h mod B. 

The idea is to count residue classes as follows, where we start with d — 0. 

Suppose d>0. For each monomial g of total degree d = deg{g), find the integer 
n such that g is new for j £ {0, . . . ,n — 1}, where we call ix'^ g new if the 
residue classes of ir^g has not already been counted. We can test the latter by 
checking whether the total degree of tt^ h mod B is at most equal to the total 
degree of h or whether n^h does not divide h. 

If we find that no monomial of total degree d contributes a new residue class, 
then we are done, otherwise we increment d by 1 and repeat this process. See 
Algorithm [T] 

Algorithm 1 Computation of length^^. Ao^E.i.v 

B = {gi{xi, . . . ,Xs^), . . . ,gr{xi, . . . ,Xs^),q} ^ Grobner basis of Id,ea,v 
m •(— ordi,(g) // q yields m distinct residue classes. 
d^O II Total degree 

T // Monomials all of whose multiples are known not to be new 
repeat 

d ^ d + 1 // Increment total degree 

V i-{g^ Y^^^ : fc, G N, X^i^i fci = d and t g for all h^T} 
1 1 No monomial of total degree d outside V can be new. 
m' <— TO 
for g do 
n ^ 

while deg(7r"(7 mod B) > d or 5 | 7r"g mod B do 

n n + 1 // TT^ij is new. 
end while 

TO <— n + TO // Get n new residue classes, 
if n = then 

T ^ T U {g} // No muhiple of g is new. 
end if 
end for 

until m = m' II No monomial of total degree d is new. 
return to 



In order to apply the results of this section we need to be able to find 

(a) an affine cover C* of C containing Dc H Ec, 

(b) ideal representations I^^i and lE,i of I?c|c» and -Eclc' • 

If the Zariski closure C of C over Spec(i?) is regular, then we can solve (a) and 
(b) by decomposing D and E into prime divisors over k. If the regular model 
has a more complicated structure, this may not be sufficient and we may have to 
decompose D and E into prime divisors over the maximal unramified extension 
fc""^ since each such prime divisor reduces to a single point on the special fiber 
of both C and C. This might be necessary because our strategy is to start with 
_R-integral ideal representations of D and E and recursively lift these through 
the blow-up process. 
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Note that we do not actually have to work over fc"''. In order to decompose 
D into prime divisors over fc"'' it suffices to consider the maximal unramified 
extension l/k contained in the smallest extension k{D)/k such that D becomes 
pointwise rational. It is possible to compute with such extensions in Magma. 

Remark 4.10. The strategy employed to decompose D and E depends on the 
curve at hand; for hyperelliptic curves there is a straightforward method to 
decompose divisors that only uses factorization of univariate polynomials as 
explained in Section O 

All of the above is trivial if D and E are pointwise fc-rational: 

Corollary 4.11. Suppose D = ni{Pi) and E — J2j '^j(Qi); where ni,mj G 
Z and all Pi and Qj are k-rational such that Dq H 0^, and Eq H contain no 
singular points of Cy . Then we have 

iy{Dc,Ec) = ^niruj min {ord„(a::i(P/) - Xi{Qj)), . . . ,ordy{xs,{Pi) - Xs,{Qj))} 
I, J 

where Pi = {xi{Pi), . . . ,Xs,{Pi)), Qj = {xi{Qj), . . . ,Xs,{Qj)) G C\ 

Proof. This follows easily from Proposition 14.81 □ 

See [20] for a similar version of Lemma [4. Ill found independently by Holmes. 

Remark 4.12. A different strategy was brought to the attention of the author by 
Florian Hess and consists in computing the intersection multiplicities of Dc and 
Ec on all affinc patches C such that Dc and Ec intersect on C*. For simplicity, 
we assume that these are and C^. We then need to subtract from the result 
all intersections that take place on both and which can be expressed as 
the length of a certain module. This approach is outlined in [IT], but it is not 
clear how it can be made practical if C 7^ C. 

4.5 Computing the correction term 



We continue to let C denote a desingularization in the strong sense of C over 
S. Suppose that the special fiber C„ is equal to X)I=o ''^i^l: where F", . . . , 
are the irreducible components of Cy. Let My denote the intersection matrix 
{iy{n,Tl,njri))^^,^^^^ ofCy. 

Suppose we are given a divisor D € Div°(C)(A:) and we want to compute a fibral 
Q-divisor $„^c(-D) = Sl=o (^i'lT-i^l having trivial intersection with Div„(C). Also 
suppose that we have found both Aly and s(Z?), where 

s{D) = {noiyiDc,T°),...,nriy{Dc,Tl)f . (6) 

We mention two possible methods here, both easily checked to be correct using 
ESI §IH.3]. 

(i) Let Af + be the Moore-Penrose pseudoinverse of My . Then we can set 

{ao,...,arf := -M+ ■ s{D). 
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(ii) (Cox-Zucker [5]) Suppose that there exists some i such that m — 1, say 
no = 1, and let be the matrix obtained by deleting the first column 
and row from My. We pick ao and 

(ai,...,a,f -.^-M'-^-s'iD), 

where s'{D) is the vector obtained by removing the first entry of s{D). 

We can now compute iy{<^yfi{D)^ Eq) easily for E e Div"(C)(A;) having support 
disjoint from D. This is simply equal to 

s{EY' ■ (ao, . . . ,ar)'^, 

where s{E) is defined as in 

It only remains to discuss how s{D) and s{E) can be computed, because can 
be computed using Magma once C has been constructed. But this is essentially 
contained in the previous section: We decompose D and E into prime divisors 
over k'^'' and then determine which components the corresponding points map to 
by lifting ideal representations recursively through the blow-up process. Since 
for any one of these divisors P all points in its support reduce to the same point 
in each step, it is easy to pick suitable affine patches covering the intermediate 
models until we find an affine patch of C containing Pq ■ The final task is the 
computation of iy[Pc,T\) for all components intersecting this affine patch 
which is easy under our assumptions. 

4.6 Computing archimedean intersection multiplicities 

In order to deal with the computation of archimedean local Neron symbols it 
suffices to consider A: = C. Let C(C) denote the Riemann surface associated 
to C. According to Section |3] we need to find an almost-Green's function with 
respect to a divisor E G Div'^(C)(C). Notice that we can write any such divisor 
in the form E — Ei — E2, where Ei and E2 are non-special, that is they 
are effective of degree g and their >C-spaces have dimension 1. It suffices to 
determine almost-Green's functions with respect to non-special divisors and 
any fixed normalized volume form on C(C). 

In order to do this it turns out to be useful to work on the analytic Jacobian, 
which we view as an abelian variety over the complex numbers. Let t be an 
element of the Siegel space f)g such that A{C) is isomorphic to the complex torus 
C9/A, where A = © rZ^. Let the map j be defined by 

j : C9 ^ Cf/A A{C). 

Moreover, we fix an Abel-Jacobi map, that is an embedding l of C(C) into 
A{C), and let Q £ Div(A) denote the theta-divisor with respect to l. Let 
S : Div(C) — ?• A denote the summation map associated to t. 
On A{C) we can find the following canonical 2-form: Let rji, . . . ,rig be an or- 
thonormal basis of the differentials of first kind on the Jacobian. Then the 
canonical 2-form is given by 

1 / - 

— (r/i A771 -f . . . + ryg Aryg) 
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and we define the canonical volume form dfi on C(C) by pulling this form back 
using t, see [HI §13.2]. The details are not important for us as the dependence 
on dfi for divisors of degree zero. 

For the next theorem, conjectured by Arakelov and proved by Hriljac, we need 
the concept of Neron functions with respect to divisors on A, for which we refer 
to Chapter 11]. We will introduce a specific Neron function in the situation 
we are interested in shortly. We use the notation Ep to denote the translation 
of a divisor E G Div(A) by a point P £ A. 

Theorem 4.13. (Hriljac) Let E e Div^(C) he non-special, let P ~ S{E) and 
E' — ([— l]*(0))p. Let Xe' be a Neron function with respect to E' . Then Xe' 
is an almost- Green's function with respect to E and dji, where d^ is the canonical 
volume form on C(C). 

Proof. See [211 Theorem 13.5.2]. □ 

The great news is that it is not difficult to find Neron functions with respect to 
0; we show below that this suffices for our purposes. 

Definition 4.14. Let 5 > 1 and a, 6 6 Q^. Let the function Oa^b on x \)g be 
given by 

Oafiiz^T)— exp f27ri f-(m + a)^r(m + a) + (to + a)"^(z + fe) 

We call 9a,b the theta function with characteristic [a]h]. 

Now let a = (1/2, . . . , 1/2), 6 = (5/2,(5 - l)/2 . . . , 1, 1/2) e Qf and consider 
da,b{z) :— 0a.biz,T) as a function on C^. 

Proposition 4.15. (Pazuki) The function Oa^ has divisor j*{Q). Moreover, 
the following function is a Neron function associated with Q: 

Xe{P) = - log \9a,biziP))\ + ^Im(z(P))^(Im(r))-i Im(z(P)), 
where j{z{P)) — P. 

Proof. This was stated without proof by Pazuki in [321 Proposition 3.2.11], but it 
is in fact rather easy to verify: It is a classical theorem (see |22l Theorem 13.4.1]) 
that the divisor of the Riemann theta function 6 =^ Sq.o is a translate by a point 
w of 8 and that 2w is the image on A of the canonical class on C. Using this it 
is not hard to see that the odd function 6a.b has divisor Then one uses 

[221 Theorem 13.1.1] to find an expression of a Neron function in terms of the 
normalized theta function 

d'a,bi^) :-0a,6(^)exp(|z^(Imr)-iz); 

the right hand side in Proposition l4.15l is equal to this expression after a straight- 
forward manipulation. 

Alternatively one can show directly that Ae satisfies the properties of a Neron 
function. □ 
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Now suppose that E = Ei — E2, where Ei,E2 € Div(C) are non-special divisors 
with disjoint support, and let Di — ^i^i{Pi) and D2 — X]f=i(Qj) t>e two 
effective divisors such that supp(i?i) n supp(Z3j) — for i,j G {1, 2}. 

Corollary 4.16. We have 
{Di- D2,Ei- E2)v 

A \ea,b{MP^)) - z(S{Ei)))6a,,{z{L{Q,)) - z{S{E2)))U 
KAMm - z{S[E2)))eaM{z[i{Q^)) - Z{S{E^))% 

d 

- 27r^ Im(z(5(£;i) - S{E2))f lTn{T)-^ Im(z(t(P,)) - z(i{Qi))), 

where for any Q E A the tuple z{Q) G satisfies i{z{Q)) ~ Q. 

Proof. Neron functions are invariant under translation of the divisor up to an ad- 
ditive constant, see [22j Theorem 11.2.1]. But according to [22j Theorem 5.5.8], 
[— 1]*(0) is just 8 translated by S'(.ft), where .ft is a canonical divisor. Hence 
the desired result follows from Theorem 14.131 and Proposition 14.151 □ 

Remark 4.17. In [2D] Holmes gives a more direct proof of Lemma [4. 161 using |53J 
§13.6/7], which relies on the theory of differentials of the third kind. 

We can use the previous result to compute intersections at archimedean places. 
In practice we need to be able to do the following: 

1) Given E G Div"(C), find non-special Ei , E2 such that E = Ei — i?2 • 

2) Compute the period matrix r. 

3) Given Pi G C(C) and r, determine 2 G such that j{z) = l{Pi)- 

4) Given t and z G C^, compute 0a.biz) — 9afi{z,T). 

5 The hyperelliptic case 



We now discuss how the methods outlined in the previous section can be com- 
bined to give a practical algorithm for the computation of canonical heights in 
the case of hyperelliptic curves. 

Suppose that C is a hyperelliptic curve of genus g defined over a field fc, given 
as the smooth projective model of an equation 

Y"^ +H{X,l)Y = F(X,l), (7) 

where F{X, Z), H{X, Z) G k[X, Z] are forms of degrees 2g + 2 and g + 1, re- 
spectively, and the discriminant of the equation ([7]) is nonzero. We will vary k 
as in general discussion of the previous sections. 

A different, but related approach to the computation of the local Neron symbols 
has been developed independently by Holmes |20j . The main difference lies in 
the computation of the non-archimedean intersection multiplicities. In |20| . 
norm maps of non-archimedean field extensions are used instead of our Grobner 
basis approach. 
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5.1 Finding suitable divisors of degree zero 

Suppose that D G Div(C)(fc) has degree zero. Then the notions introduced 
in Section 14.11 are all well-known: The reduction process is part of Cantor's 
algorithm for the addition of divisor classes introduced in ^ ; here the divisor A 
used for reduction is equal to (oo) when we have a A:-rational Weierstrass point 
oo at infinity and is equal to (ooi) + (002) when there are two branches ooi, cxi2 
over the singular point at infinity in the projective closure of equation ([7]). 

In the former case Lemma l4.1l savs that the reduction process yields the unique 
effective divisor D such that 

D ^ D + r{oo), 

where < — r = degD < g and deg(D) is minimal. In the latter case it turns 
out that when g is even we can still find a unique D of minimal nonnegative 
even degree —r<g such that 

if we impose further conditions on its ideal representation. Conversely, if g is 
odd we might have to take reductions of degree 5+1 into account and these are 
not unique. However, uniqueness of the reduction is not an essential property 
in our applications and so we shall not discuss it any further. 

A possible ideal representation of a reduced effective divisor D is given by the 
Mumford representation which we now recall briefiy. 

If we view C as embedded in weighted projective space with weights 1, <? + 1, 1 
assigned to the variables X,Y, Z, then it is given by the equation 

y2 + H{X, Z)Y = F{X, Z). 

An effective divisor D of degree d < g + 1 corresponds to a pair of homogeneous 
forms {A{X, Z),B{X, Z)), where A(X, Z) and B{X, Z) have degrees d and g + 1 
respectively, such that D is defined by 

A{X, Z) = = Y -B{X,Z) 

and we impose the additional condition that 

A{X, Z) I B(X, Zf + H{X, Z)B{X, Z) - F{X, Z). 

First suppose that there is a unique Weierstrass point 00 at infinity in C{k). 
Then any nonzero effective divisor D = X]j=i(^i) that is reduced along (00) 
has degree d < g and cannot contain 00 in its support. Hence we can safely 
dehomogenize in order to represent D and so we may take 

Id = {a{x),y- b{x)), 

where a{x) = A{x,l) and b{x) — B{x,l), for its ideal representation. More 
concretely, we have 

n 

a{x) = l[{x-x{P,)) 
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and b(x) has minimal degree such that 

b{x{P,)) ^ y{P,) forj = l,...,d. 



Conversely, suppose that there are two points cxdi,oo2 at infinity. Suppose 
that D is reduced along (ooi) + (002). If supp(£') does not contain a point at 
infinity, then we can dehomogenize as before to find an affine representation. 
If this does not hold, say ooi € supp(D), then necessarily 001,002 6 C{k) and 
002 ^ supp(I?). This case is more subtle, because we cannot tell the multiplicity 
of ooi in Z? from its dehomogenized form. For our applications it suffices to treat 
the affine and the infinite part of D separately. Hence this complication does 
not cause any trouble. 

Now let fc be a global field, and let the divisor Doo be defined by 2(oo) if there 
is a unique fc-rational point at infinity and by (ooi) + (002) otherwise. Also 
suppose d is even and 

where D = reduced along Doo such that no Pi is a point at infinity 

or a Weierstrass point. Then we can always use ni = 1 and n2 = — 1 in the 
method introduced in Section HTTl this is due to Holmes, see [20]. Namely, if we 
apply the hyperclliptic involution 

to the points Pi, then we have 

i=l 

If we move this by the divisor of a function x ~ C, where C g fc is such that 
x{Pi) ^ ( for all Pi, then we find 

supp(D) n supp(£:) = 0, 

where E = D' + d/2div{x — (). This corresponds to choosing A — Doo and 
A' = D{() in the method outlined above, where D{C,) = div(x — C) + ^00 • 

Instead of computing {D,D), we can now compute 

h{P) = -{D,D) ^ {D,-D) ^ {D,E). 

If we have 

D = D-^Doo, 

where 

d' 

D = ^{Pi) + UooioOl) 
i=l 

is reduced along D^o — (ooi) + (002), such that d — d' + Uoo and all Pj are 
affine non- Weierstrass points (see Section [5TT|) . then we also have to move D 
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away from ooi using a function x — C, where x{Pi) ^ C' 7^ C foi' alH = 1, . . . , d'. 
The computation becomes 

d' , d' , 

-{D,D) = {J2{P^)+r^oo{^l) ~ -D{C).Y.^P- ) ^ n^{^^) - -D{0) 

i=l i=l 

and poses no additional problems due to the bilinearity of the local Neron sym- 
bol. 

What if there is a unique rational Weierstrass point oo at infinity and d is odd? 
In that case we use 

d 

D' = 2Y,{P~)~dD^^ -2D 

i=l 

and compute 

h{P) = -{D, D) = {D, -D) = ]^{D, E), 

where E = D' + ddiv{x — and C is as above. Note that we can still use the 
reduced Mumford representation, because we have 

d 

{D,E) = 2{D,Y,{Pr))-d{D,D{0). 

i=l 

Finally, if supp(D) contains an affine Weierstrass point, then we simply compute 
h{P) = -^h{nP) such that nP has a reduced representation not containing an 
affine Weierstrass point. 

5.2 Determining relevant non-archimedean places 



Suppose that A: is a global field. Our curve C is covered by two affine patches 
and C^, where 

■.y^ + Hix,l)y^Fix,l) (8) 

and 

: + Hil,z)w ^ F{l,z). (9) 

It follows from the discussion at the end of Section 14.21 that we can assume 
Ok to be Euclidean. Suppose that D and E are Id.i = {a{x),cy — b{x)) and 
Ie,i = {a'{x), c'y — b'{x)) on , respectively (where we have multiplied all poly- 
nomials by the least common multiple of the denominators of their coefficients, 
if necessary) . Then we need to compute a Grobner basis of 

Id,es = {y^ + H{x, l)y - F{x, 1), a{x), a'{x),cy - b{x),c'y - b'{x)) 

and factor the unique element qD,E,i of Ok appearing in this basis. 

Now suppose that v G M° satisfies iv{Dv,CyEv,c) > 0, where C -> C" is a 
desingularization in the strong sense and that the points of intersection do not 
map to the closure of in C. Any such v must satisfy v{ad) > and v{a'^,) > 0, 
where ad and a'^, are the leading coefficients of a{x) and a'(x), respectively. So 
instead of computing a Grobner basis of Id,e,2, we can factor gcd(ad, a'^) which 
is usually much easier than factoring {qD,E,2)- This simplification can make a 
big difference in practice. 



19 



5.3 Computing non-archimedean intersection multiplici- 
ties and the correction term 



Let k denote a non-archimedean local field and let tt be a uniformizing element. 
Let D £ Div(C)(fc) be effective such that an ideal representation of D on is 

Id,i = {a{x),y - b{x)), 

where a{x), b{x) G k[x] and we have deg(a) < g and deg(6) < 17 + 1 as in Section 

In order to use Proposition 14.81 to compute non-archimedean intersection mul- 
tiplicities, we need to be able to decompose divisors into prime divisors over 
unramified extensions. The main point distinguishing the hyperelliptic case 
from the general situation is that we can decompose divisors by factoring uni- 
variate polynomials over non-archimedean local fields; we show how this can be 
used in this section. 

Note that factoring univariate polynomials over p-adic fields and Laurent series 
over finite fields is implemented in Magma (following work of Pauli [311). 

We first deal with the case C — C and use the affine cover C = U , where 
is the Zariski closure over S of the affine curve C* defined as in Section [5T2] 

We can factor a{x) = ai{x)a2{x), where 02(2;) is constant modulo tt and ai{x) £ 
R[x]. This corresponds to a decomposition D = Di + D2, where -Di.c lies in 
and D2fi lies in C^. More precisely, we have 

IdiA = (ai(x),y - 61 (x)), 

where bi{x) = b(x) mod ai{x). In order to use Proposition l4.81 we need 61 (x) £ 
R[x]. Suppose that ai{x) is irreducible (otherwise factor ai{x) into irreducibles) 
and that bi{x) ^ R[x]. If Di e does not have a singular point of the special fiber 
Cy in its support (for instance, if ai is unramified), then we can safely extend k 
by a root of ai and work over this extension. Repeating this process, if necessary, 
leads to a finite extension k' of k such that Di e splits into prime divisors over 
k' that have _R'-rational ideal representations, where R' is the ring of integers 
of k'. 

Now suppose that ai{x) reduces to (x — a)™ mod tt, where a is the x-coordinate 
of a singular point of Cy and m > 1. In general we cannot extend the field 
by a root of oi, because there may be points in the support of Di^c that are 
not regular over this extension. But because of the special shape of ai, we 
can simply use the i?-rational ideal representation (ai(x), tt'^j/ — b[{x)), where 
bi{x) = TT~''b[{x), and bi{x) € R[x] has a unit among its coefficients. Note that 
this approach does not always work for more general ai{x). 

If we have C =^ C, then we simply start by factoring a{x) into irreducibles over 
fc""^. Assuming that ai{x) is one of the irreducible factors, we lift the ideal 
representation Idi,c = (ai(x), Tr^y — b'i{x)) recursively through suitable blow- 
ups until we arrive at a suitable affine patch where the intersection multiplicities 
can be computed using Proposition 14.81 As explained in Subsection 14. 5[ this is 
also sufficient to compute the correction term. 
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5.4 Computing archimedean intersection multiplicities 



In order to compute archimedean intersection multiplicities, we need algorithms 
for steps l)-4) introduced at the end of Section WM 

For hyperelliptic curves, steps 2), 3) and 4) have all been implemented in Magma 
by van Wamelen. An earlier version of the implementation using Mathematica 
can be found on van Wamelen's homepage [4^. The routines there only work 
for genus 2 curves, but the general algorithms in Magma work similarly. 

It is important to note that step 4), the computation of 0a,b, is done via ap- 
proximation using the definition, in particular it can be used to compute 9a,b 
provably up to desired precision. 

We discuss step 1). Here we want to find, given P, Q G v4, divisors Di, D2, Ei 
and E2 such that 

(a) [Di -D2] = P and [^i - E2] = Q, 

(b) Di, D2, El, E2 are effective and have pairwise disjoint support, 

(c) El and E2 are non-special. 

We can allow ourselves more freedom, and only require that (a) holds for some 
multiple nQ, due to the bilinearity of the local Neron symbol. For simplicity we 
only discuss the case of a unique point at infinity, the other case being similar 
with a few minor subtleties if g is odd. We pick Di :— D and D2 := d{oo) if 
P is represented hy D — d{oo) and D has affine support. Suppose that nQ is 
represented by En — 5(00), where En is non-special and has affine support such 
that En and E!^ have support disjoint from D, where E'n is the result of the 
hyperelliptic involution applied to the points in the support of En ■ Then 2nQ 
is represented by En ~ En and we choose Ei := £"„ and E2 := E!^. . 

With these choices, at most d + g applications of the Abel-Jacobi map and at 
most 2d applications of the theta-function 0a,b are required in order to compute 
{Di — D2,Ei — E2)v for an archimedean place v, essentially because we have 
t(oo) 0. 

Now let C G /c* be as in Section [01 We are actually interested in computing 

{D-d{^),E-e/2D{C)),, (10) 

so we compute a function /3 e k{C)* such that 

div(/3) =Ei-E2- 2nE + ndD{C) 

See [19] for an algorithm that computes /3. Using properties (a) and (c) of 
Proposition 12 . 81 we can compute (|T0)) . 

Notice that in contrast to the non-archimedean case the running times of steps 
3) and 4) do not crucially depend on the heights of the points in the supports 
of the respective divisors, since we work with the complex uniformization. 
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6 Examples 



In this section we provide a hyperelliptic example of a regulator that was 
computed using the algorithm outlined in the previous sections. Moreover, we 
shall discuss, at least in the case of hyperelliptic curves, how the running time 
changes as we increase 

(a) the genus of the curve; 

(b) the size of the coefficients of the point. 

All times are in seconds, unless noted otherwise. For the computations we have 
used 50 digits of p-adic and 30 digits of real and complex precision. 

We first use the Magma-implementation of our algorithm to compute the regula- 
tor of the Jacobian of a hyperelliptic genus 3 curve up to an integral square. 

Example 6.1. Let C be given by the smooth projective model of the equation 

= X{X - 1){X ~2){X - 3){X -6){X- 8){X + 8). 

The curve C is a hyperelliptic curve of genus 3, defined over Q. A quick search 
reveals the following rational non-Weierstrass points on C. 

(-2, ±240), (4, ±48), (-6, ±1008) 

Let A denote the Jacobian of C; obviously its entire 2-torsion subgroup is defined 
over Q. In order to bound the Mordell-Weil rank of A we compute the dimension 
of the 2-Selmer group of A over Q using Magma. This dimension is equal to 3 
and hence we get an upper bound of 3 on the rank. If Pi, . . . , P„ G A, then we 
denote the regulator of Pi , . . . , P„ by 

Reg(Pi, . . . , P„) = det ((P„ P,)nt),,, • 

We want to compute the regulator Reg(P, Q, R) of the subgroup G of A{Q) 
generated by the points 

P = (-2, -240) - (oo) 
Q ^ (4,-48)-M 
R = (-6, 1008) - (cx)). 

The discriminant of C factors as 2''^"3^^5^7"'ll^. We first find regular models 
at the bad primes 2, 3, 5, 7 and 11. All computations in this example were 
done using Magma on a 1.73 GHz Pentium processor. It turns out that all 
computed regular models are already minimal; we list the number of components 
of the special fiber of the respective regular model, the (geometric) group of 
components of the Neron model and the time it took to compute the regular 
model in Tabled] 

After this preparatory step we now compute the entries of the height pairing 
matrix. The results and timings can be found in Table [H 
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prime 


^ of comps. 




time 


2 


14 


(Z/2Z)^ 


1.95 


3 


9 


(Z/2Z)3 X Z/4Z 


0.35 


5 


4 


(Z/2Z)3 


0.23 


7 


3 


(Z/2Z)2 


0.29 


11 


2 


Z/2Z 


0.10 



Table 1: Regular model data 



S e A{Q) 


h{S) 


time 




P 




1.90008707521104082692048090266 


23.10 




Q 




1.15261793630905629106514447088 


19.76 




R 




2.90090831616336727010940214290 


20.96 


p 


+ 


Q 


2.36481584203715381857836835238 


19.95 


p 


+ 


R 


5.51584078564985349844572029952 


20.67 


Q 


+ 


R 


5.74901893484137170755580219303 


21.22 



Table 2: Canonical height computations 



Using these results, we find 

Reg(P, g, R) = 4.28880986177463283058861934366. 

We can test our findings by computing Reg(nP, rnQ, IR) for several integral 
values of n, m, /. In all cases we get the relation 

Reg(nF, mQ, lR)/Rcg{P, Q, R) = n^m^f 

up to an error of less than 10~^^. As Reg(P, Q, R) is clearly non-zero, we know 
that G is a subgroup of finite index and hence Reg{P,Q,R) equals Reg{A/Q) 
up to a rational square. 

Next we want to illustrate the behavior of the running time of our algorithm. We 
have refrained from a formal complexity analysis, mostly because the algorithm 
uses several external subroutines, such as the computation of regular models 
and of theta functions, whose complexities have not yet been analyzed. 

In the case of zero-dimensional ideals of polynomial rings over fields, the com- 
plexity of a Grobner basis computation can be shown to be polynomial in £)", 
where D is the maximal degree of the elements of the basis we start with and 
n is the number of variables. See for a summary of results regarding com- 
plexity of Grobner basis computations. In particular this holds for Faugere's 
F4-algorithm [Tl], used for instance by Magma (over fields and Euclidean rings). 
This result can be extended easily to the case of polynomial rings over Euclidean 
domains, provided we have fast algorithms available for the linear algebra com- 
putations in the i<"4-algorithm, such as those implemented in Magma. So the 
Grobner basis computations do not cause any trouble in practice, since the way 
regular models are computed in Magma ensures that the number of variables does 
not grow. 

Indeed, the running time of the algorithm is usually dominated by the various 
analytic computations required for the archimedean local Neron symbols. They 
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depend exponentially on the genus; the curve of largest genus that we have 
been able to compute with has genus 10, see Example 16.21 below. If the genus 
is not very large, but the size of the coefficients of the point P £ A{k) that 
we want to compute the canonical height of is, then it turns out that the main 
bottlenecks are usually the factorizations alluded to in Section l42l recall that 
these are required in order to find out which places can lead to non-trivial non- 
archimedean local Neron symbols. See Example 16.31 



d 


genus 


HP) 


act 


nact 


5 


2 


1.20910894883943045491548486513 


3.51 


0.33 


7 


3 


1.31935353209873515158774224282 


6.70 


0.34 


9 


4 


1.39237255678179422540594853290 


12.65 


0.87 


11 


5 


1.44187308116714103129667604112 


32.30 


1.67 


13 


6 


1.47679608841931245229396457463 


120.51 


2.99 


15 


7 


1.50265701979128671544005708236 


791.14 


5.17 


17 


8 


1.52254076352483838532148827258 


4729.03 


8.95 


19 


9 


1.53829882683402848666502818888 


62535.55 


14.20 


21 


10 


1.55109127084768378637549292754 


280731.59 


21.35 



Tabic 3: Canonical heights in a family 



n 


h{nP) 


act 


nact 


factt 


digits 


1 


0.19809838973401855248161508134 


2.35 


0.02 


0.00 


1 


2 


0.79239355893607420992646032538 


2.38 


0.02 


0.00 


1 


3 


1.78288550760616697233453573211 


4.06 


0.13 


0.00 


1 


4 


3.16957423574429683970584130153 


3.36 


0.11 


0.00 


1 


5 


4.95245974335046381204037703364 


2.91 


0.10 


0.00 


1 


6 


7.13154203042466788933814292846 


3.39 


0.08 


0.00 


3 


7 


9.70682109696690907159913898594 


3.40 


0.06 


0.00 


6 


8 


12.6782969429771873588233652061 


3.36 


0.34 


0.05 


9 


9 


16.0459695684555027510108215890 


3.31 


0.29 


0.01 


11 


10 


19.8098389734018552481615081345 


3.41 


0.95 


0.64 


16 


11 


23.9699051578162448502754248428 


3.33 


0.37 


0.08 


19 


12 


28.5261681216986715573525717137 


3.45 


0.47 


0.11 


21 


13 


33.4786278650491353693929487474 


3.32 


0.34 


0.09 


21 


14 


38.8272843878676362863965559437 


3.42 


196.87 


196.50 


30 


15 


44.5721376901541743083633933028 


3.37 


0.53 


0.20 


38 


16 


50.7131877719087494352934608245 


3.39 


0.90 


0.25 


42 



Table 4: Canonical heights for multiples of a point 



All computations for the following two examples were done using a 3.00 GHz 
Xeon processor. 

Example 6.2. Consider the family 

Cd:y'^^x'^ + 3a;2 + 1 

for d e {5,7,9,11,13,15,17,19,21} and let P = [(0,1) - (0,-1)] e AdiQ), 
where Ad is the Jacobian of Cd- We compute h{P) and record the running time 
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for both the archimedean and the non-archimedean computations. See Table 
[31 where nact and act denote non-archimedean and archimedean computation 
time, respectively. This example illustrates the exponential dependency on the 
genus. 

Example 6.3. Next we consider C : y"^ — — a:'^ + 1 and let P G A{Q) = 
[(0,1) — 0O+] G Jac(C)(Q), where oo+ is the point at infinity such that the 
function y/x^ has value 1 at 00+. The curve C has bad reduction at 2. We 
use Magma to compute a regular model at 2; this takes 1.97 seconds and yields 9 
irreducible components. 

We compute the canonical heights of positive multiples of P £ A{Q) and record 
running times. The results are in Table |4] and we see that we have h{nP) — 
n?h{P) for all n G {1, ... , 16}. Here nact and act have the same meaning as in 
Table [3j factt denotes the time spent on integer factorization and digits denotes 
the number of digits of the maximal height of the coefficients of the polynomials 
in the Mumford representation of nP. 

For n > 17, the time spent on factoring increases dramatically. For instance, 
the factorization needed for n = 18 takes about 62 hours. The largest multiple 
for which we are able to compute h{nP) is n = 21, where digits= 79. 

7 Outlook 



It is now possible, using the Magma-implementation of the algorithm described 
in this work (see [22]), to compute canonical heights on Jacobians of hyper- 
elliptic curves defined over number fields. There is work in progress on most 
of the applications outlined in the introduction. Some can now be tackled in 
a straightforward way, such as the computation of regulators up to integral 
squares, which can be used to gather numerical evidence for the conjecture of 
Birch and Swinnerton-Dyer as in [16j . some require more work to be done first, 
such as the computation of generators of the Mordell-Weil group. An algorithm 
for the latter is presented in [37!, but in order to apply it, one also needs a suit- 
able naive height combining the properties that we can list all points of naive 
height up to some bound and that the difference between the two heights can 
be bounded effectively. Holmes has recently come up with a good candidate for 
such a naive height, the details will appear in his upcoming PhD thesis at the 
University of Warwick. For the genus three case, there is recent work of StoU 
[39] solving this problem. 

We now sketch some possible directions for further research regarding the canon- 
ical height algorithm itself. First, our algorithm works for any global field and 
hence it should not be too difficult to implement it for hyperelliptic curves de- 
fined over global function fields. In fact, some problems disappear because of 
the absence of archimedean places. More importantly, it would be interesting 
to extend our algorithm to the case of non- hyperelliptic curves. Here, there are 
essentially two problems: 

(i) How can we decompose divisors into prime divisors? (see Section [4.21) 
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(ii) How can we implement the analytic steps 2) - 4) introduced at the end of 
Section mF 

There are 3 approaches to problem (i). If we could factor multivariate polyno- 
mials over non-archimedean local fields, then (i) would be solved, but such an 
algorithm has not been implemented to the author's knowledge. In favorable 
situations it may be possible to use ideal representations similar to the Mumford 
representation of hyperelliptic curves and thus decompose divisors using factor- 
ization of univariate polynomials. An ideal representation resembling Mumford 
representation has been proposed in |34] for smooth plane quartics. Finally, 
it might not be necessary to decompose divisors at all, if we could make the 
approach mentioned in Remark l4.12l and described in |41) work in our situation. 

Of course, problem (i) disappears whenever we deal with divisors having point- 
wise /c^-rational support for each relevant non-archimedean local field ky. We 
have used this to compute all non-archimedean local Neron symbols necessary 
for the computation of the regulator (up to an integral square) of the Jacobian 
of a non-hyperelliptic curve of genus 4 without special properties, see [281 Chap- 
ter 6]. This curve plays a major part in [38) . where it is shown, assuming the 
first part of the conjecture of Birch and Swinnerton-Dyer, that there are no ra- 
tional cycles of length 6, an important result in arithmetic dynamics. Our goal 
was to verify the second part of the conjecture of Birch and Swinnerton-Dyer 
up to an integral square for this particular curve, a challenge problem posed by 
Stoll in [3H]. 

Regarding problem (ii), an extension of van Wamelen's algorithms to the non- 
hyperelliptic case would suffice. This does not appear to be particularly difficult, 
but we have not attempted it. The main obstacle is the choice of a basis of 
holomorphic differentials, which will probably depend on the class of curves 
under consideration (for hyperelliptic curves there is, of course, a canonical 
choice and this is used by van Wamelen's algorithms). 

If we are willing to allow non-rigorous numerical integration methods, then all 
of the relevant algorithms have been developed (see [10j[5l[Tl]) by Deconinck et 
al. for general compact Riemann surfaces. However, these algorithms are not 
suitable in order to actually prove that we have computed the correct canonical 
height up to given precision. 

In any case, Deconinck and his collaborators implemented their algorithms in 
Maple in a package called algcurves. Unfortunately, the Maple developers have 
since decided to change some of the functions that algcurves uses, in the process 
destroying some of the package's crucial functionality. Deconinck's group are 
currently working on a long-term project to rewrite all necessary routines in 
Sage. 

Finally, it would be interesting to formally analyze the complexity of our algo- 
rithm. While knowing that the dependence on the genus is exponential due to 
the necessity of computing theta-functions, we first need to analyze the com- 
plexity of Magma's desingularization algorithm and the analytic algorithms that 
we use before more can be said. 

Remark 7.1. Let A be the Jacobian of a smooth projective curve C defined over 
a number field and let p be a prime such that A has good ordinary reduction at 
all V I p. Coleman and Gross [S] have constructed a pairing taking values in Qp 
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between divisors D,E on C of degree which can be decomposed into a sum of 
local height pairings hy{D, E) over all non-archimedean v G Mk- 

They show that this pairing respects linear equivalence and coincides with the 
p-adic height pairings on A constructed by Schneider [3S] and Mazur-Tate p7] . 

If u I p, then hy{D,E) can be expressed in terms of Coleman integration. An 
algorithm for the computation of hy{D, E) was introduced by Balakrishnan and 
Besser in g] , see also Chapter 8] . Computing the local height pairings at v ] p 
is equivalent to computing the local Neron symbol at v (cf. [21 §2]). Hence we 
can combine the results presented in this work with the method of Balakrishnan 
and Besser, leading to the first algorithm to compute p-adic heights for g > 2. 

One can define p-adic regulators similar to the classical case (see Section [S]) . 
Together with Balakrishnan we have computed the p-adic regulator for all but 
one of the modular abelian surfaces considered in [TB] for all good ordinary 
p < 100. This gives rise to an extension of the p-adic Birch and Swinnerton- 
Dyer conjecture for elliptic curves due to Mazur, Tate and Teitelbaum to the 
case of modular abelian surfaces (cf. [S] Conjecture 9.1.4]). 
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